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We calculate the tt, p, uj, and 7 coupling constants between the heavy meson doublets (1~,2~) 
and (0~, l~)/(0''', 1^) within the framework of the light-cone QCD sum rule at the leading order 
of heavy quark effective theory. Most of the sum rules are stable with the variations of the Borel 
. parameter and the continuum threshold. Then we calculate the strong and electromagnetic decay 

' widths of the (1~,2~) D-wave heavy mesons. Their total widths are around several tens of MeV, 

which is helpful in the future experimental search. 
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I f'. I. INTRODUCTION 

^ : 

^p^. Heavy quark effective theory (HQET) [ij is a framework which is widely used to study the spectra and transition 
O i' amplitudes of heavy hadrons containing one heavy quark. In HQET, the expansion is performed in terms of I/ttiq, 
1^ where mg is the mass of the heavy quark involved. At the leading order of 1 /toq , the HQET Lagrangian respects the 
heavy quark flavor-spin symmetry, therefore heavy hadrons form a series of degenerate doublets. The two members 
' in a doublet share the same quantum number j;, the angular momentum of the light components. The ji = i 
S-wa.ve doublet (0^,1^) is conventionally denoted as H and the ji = P-wave doublets (0+, l+)/(l+, 2+) are 
, conventionally denoted as S/T. We denote the ji = |/| D-wave doubtlets 2~)/(2~, 3~) as M/N. 

Shifman-Vainshtein-Zakharov (SVZ) sum rules is a nonperturbative approach used to determine hadronic pa- 
rameters such as the hadron mass. The vacuum expectation value of the T product of two interpolating currents is 
considered in this approach. After performing the operator product expansion (OPE), one obtains sum rules which 
relate the hadronic parameters to expressions containing vacuum condensates parameterizing the QCD nonperturba- 
tive effect. In the late 1980s, light-cone QCD sum rules (LCQSR) Q was developed to calculate various hadronic 
^) ] transition form factors. Now the OPE of the T product of two interpolating currents sandwiched between the vacuum 
T-H . and an hadronic state is performed near the light-cone rather than at a small distance as in the conventional SVZ 

■ sum rules. 

, The p coupling constants gB*Bp and go* Dp were calculated with LCQSR in full QCD in Ref. The couplings 

■ gH*H*p, fH*H*p, gHHp, and fH*Hp were calculated in full QCD in Ref. @. Their values in the limit toq — 00 are also 
discussed in this paper. The p coupling constants between the three doublets H, S, T and within the two doublets 

- - I H, S are systematically studied with LCQSR at the leading order of HQET in Ref. The tt coupling constants 
between the S'-wave and P-wave heavy mesons have been studied using QCD sum rules or/and LCQSR in Ref. Q- 
The TT coupling constants between M/N and H/S/T are calculated with LCQSR at the leading order of HQET in 
Ref. [8, 9]. The radiative decay between H, S, and T are studied using the light-cone QCD sum rule at the leading 
order of HQET in Ref. 10]. In Ref. the radiative decays of Dlj{2317) and Dsj{2460) are studied using LCQSR 
approach. 

In this work, we use LCQSR to calculate the tt, p, uj, and 7 coupling constants between the doublets M and H/S. 
Because of the covariant derivative in the interpolating currents of the M doublet, the contribution from the 3-particle 
light-cone distribution amplitudes of tt, p, w, and 7 have to be included. We work in HQET to differentiate the two 
states with the same value and yet quite different decay widths. The interpolating currents Jj'pj"' adopted in 
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our work have been properly constructed in Ref. [l^. They satisfy 

(0|j;:-J^(0)|j',P',jn = /P,,%'<5PP'<5,,,.,7— ^ (1) 

in the hmit mg — oo. Here 7y"i '"j is the polarization tensor for the spin j state, v is the velocity of the heavy 
quark, g"^ = — v"v^ , S denotes symmetrizing the indices and subtracting the trace terms separately in the sets 
{ai---aj) and 

II. SUM RULES FOR THE n COUPLING CONSTANTS 

We shall perform the calculation at the leading order of HQET. According to Ref. [Ill, the interpolating currents 
for the doublets (0^,1"), (0+,l+), and (1~,2-) read as 



= y^^.(-»)75 i^ir-D?' + ir-Dr - ^^r^"^ a} ? , (3) 

where hy is the heavy quark field in HQET, 7^ = 7'* — ■Ou'', 2?^ = 2?'^ — (P • v)v^^ , = 5^"^ — u^u'^, and u'^ is the 
velocity of the heavy quark. 

We consider the tt decay of M2 to Hi to illustrate our calculation. Here the subscript of M{H) indicates the spin 
of the meson involved. Owing to the conservation of the angular momentum of the light components in the limit 
iTiQ — > 00, there is only one independent tt coupling constant between doublets M and H. We denote it as g^j^ju^ 
where p and the number following it indicate the orbital and total angular momentum {l,jh) of the final tt meson 
respectively. g^]ij^ can be defined in terms of the decay amplitude A^(M2 Hi + n) as 

MiM2 -> i?i + tt) = [er^qr - Igr'^' {e* ■ qt)]gfj,H,. , (4) 

where rj and e* denote the polarization tensors of the initial and final heavy mesons respectively, q is the momentum 
of the TT meson. The transversal tensor are defined as = — {e ■ v)v^, q^ = q'^ — {q ■ v)v^, and g^'^ = g^'^ — v^v'^ . 
/ = 1, 1/a/2 for the charged and neutral tt meson, respectively. 

To obtain the sum rules for the coupling constants g^f^ffiTr' consider the correlation function 



d^xe-^'^-- (7r(g)|r{ Jf , {Q) jI'^-J {x)m 



I 



\{gT'qT^9T'qT)-\9T'''q^ 



G1U,.(-:-'), (5) 



where a; = 2u • fc, cj' = 2u • (fc — q). At the leading order of HQET, the heavy quark propagator reads as 

(0|r{/z„(0)/i„(x)}|0) = / dt6\-x - vt) . (6) 
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The correlation function can now be expressed as 



75 



■;9t 



{7T{q)\q{x)m\0) 



(7) 



It can be further calculated using the light-cone wave functions of the tt meson. To our approximation, we need the 
2- and 3-particle light-cone wave functions. Their definitions are collected in Appendix \X\ 
At the hadron level, G^J^^^^(aj, oj') in ^ has the following pole terms 



/-4/2/-,3/25M2//i7r 



(2A_,i/2-Lj')(2A_,3/2-cj) 2A_^i/2-w' 2A_,3/2 



(8) 



where A_ x/2 = — ^nq, A__3/2 = "rnM — mq, /_^i/2, etc. are the overlap amplitudes of their interpolating currents 
with the heavy mesons. 

G^j^j^fj_^p{ui,ui') can now be expressed by the tt meson light-cone wave functions. After the Wick rotation and the 
double Borel transformation with oj and oj', the single-pole terms in ([8]) are eliminated. We arrive at 



^-,3/2+^-. 1/2 



48' 



-/J -12 K(So) - (u0.)'(uo)] T'hi^) 



ruu + rud 



T 



'(uo) - (uAYiuo) + 8B[2l(So) - 8B[il(i2o) -I- 8(mB)W(?2o) - 16Vf'°\uo) - 16yt[f''"'(Mo) 



,[0,0] 



(9) 



where fn{x) — 1 — e ^X^I-^o^*^/^' ^® continuum subtraction factor, and Wc is the continuum threshold, uq = 
Ti/{Ti + T2), T = TiT2/{Ti + T2), and uq = 1 - uq. Ti and T2 are the two Borel parameters. We have employed 
the Borel transformation B'^e°"^ — 6{a — 1/T) to obtain In the above expressions, we have used the functions 
J"['^l(uo) and J"["'''l(uo) which are defined in Appendix [Bl 

The TT coupling constant between doublets M and S/T/M can be defined similarly: 



Tin fPaiqv 012 d2 



M{M2 ^Si+tt) 
M{Ah -> Ti + ^) - 7(77 • e;)gt>^j,^^ + I 



9? 



d2 
9MiTi-k 



M{M2^Ah+TT) = 2/77^,, 



-,9t (e -It) 



+/ry„,„, ^ q^qr {e* • g,) " | P^r^C^ + 9?'"' {e* ■ qt)] \ .g^^.M^. ■ 



pi 

9m2Mitt 



J3 



(10) 



Here the vector notations in the Levi-Civita tensor come from an index contraction between the Levi-Civita tensor 



and the vectors, for example, e 
the coupling constants in Eq. (ITUl) are 



q^vg. The Levi-Civita tensor is defined as e 



0123 



1. The sum rules for 



-,3/21^^^+, 1/2 



24' 



~3m. 



T ' mu + ma 
(uo) - (uA)(uo) + 16yt[|"''°i(uo) + IQAZ^'°\uo) 



,3/2+^^+, 3/2 



^0^2) ^ (^(^ _ z.)0p)(2)(So) + (1 - a2r)[3'"l("o) 



-(a3r)[^'il("o)] TV2(^)-^' 



^A'(72o) + -^(^(1 - u)A)(3)(So) + bW(uo) + 
8 16 



2»n"oj 



- u)B)'(uo) + {u{l - u)ct>„y{uo) + (2^11 + 2A± + Vy + V±)^^'^\uo) + ((1 - a2){A\\ + V±)f'^\uo) 
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- u)0p)(mo) + 



ml - m. 



ud 



6ml 



16' 



+B[3l(,io) + (^11 - V|| - Vi)[-i'°](«o) - ((1 - a2){Ai\+V^)f'"Kuo) + (a3(^|| + Vi))[°'il K) 



9MiTi7Tf+,if-,i<^~ 



,3/2 + '^+, 3/2 



^(«(l-t.),/..)'MTVi(S)- 



T 4TO„d 



- u)(j)p){uo) + 
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ml — m 



6ml 



^0aM + ((i-a2)r)[i'°iM 



-(a3r)^iJMjT/o(^^)+ 2 



- «)A)'(iio) - l{u{l - ti)B)W(fZo) +B[3l(«o) + (Oi - ^)B)Pl(wo) 



+(^|| +2V|| +2Vi)[-i'°lK) - ((1 - a2)M|| + Vi))[°'01(«o) + (asC^n + V±))[°'°1K) 



pi /•2 

jT^m^ud 



,3/2 



80\/6 



T 240\/6m„d 



(«(1 - «)(/),)(2)(uo) + 6((1 - a2)r)[''°l(uo) 



20V6 



|a(So) - - «)A)('HSo) - (^^(1 - u)^^){uo) 
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-4(^11 +^^)["^ol(^io) + ((1 - a2)(^|| +^x))l''°l(«o) - (a3(^|| - 5(V|| + Vi)["^°lK) 



Tfo{—) + ■^^"^^ 
T 60\/6m„<j 



mi. — m: 



- «)<^<.)(«o) - 6((1 - a2)r)[°'°l(uo) + 6(a3r)[°'^l(«o) 



Am^ 
5^6 



^(u(l - u)A){uo) - (^11 + ((1 - a2)(^|| + - (a3M|| +^±))'-'''l(«o) 



1 



,3/2 



/3 j-2 
S'MgMiTr-'-,!^ 



2mw 



m' - m: 



6m?. 



HI (So) - ((1 - a2)r)[°'°iK) 



- ti)A)(wo) - (^11 + ((1 - a2)M|| +^±))'-''°l(«o) 



+(a3r)[°'^]K)J +\/6^m2 

-(a3(^ll+^i))[-i'ilM" 
with m„d = m„ + m^. The tt coupUng constants of the other decay channels are defined as 
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(11) 



M{Mi^ Hi +7r) = /ie''^*'''^ 



M{Mi Si+n) 

M{M2 So + tt) 

M{Mi -^T2+n) 
M{M2 -s- Ti +7r) 
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9M2T2TT 
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-,9t 9t 



9t 9t 



d2 

9M2T2TT 



M(Mi Mi+tt) 



M{M2^Mi+n) = 2I'na,a2 



pi 

9M2MiTT 



MiAh^M2 + 7r) = ^nVo.,a2Vkf>2^"'^""'9r^'g'M2M2 



a2^02 It 0202 

It It "g"yt 



/3 

9m2M2TT' 



(12) 



We notice that 
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/3 
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(13) 

These relations are consistent with the HQET expectation at the leading order. 

To get the numerical values of the six independent coupling constants, we need the following mass parameters A's 
and /'s, the overlap amplitudes of these interpolating currents: 

(0-,l-)fl3] {0+,l+)\U\ {l+,2+)\Ml (l-,2-)[9] 

A 0.5 GeV 0.85 GeV 0.95 GeV 1.42 GeV 

/ (0.25 ± 0.04) GeV^/^ (0.36 ± 0.10) GeV^^^ (0.26 ± 0.06) GeV^/^ (0.39 ± 0.03) GeV^^^ 

The TT decay constant = 131 MeV. /x^ = ml/{mu + md) is given in Ref. [ll]: A^tt = (1.573 ± 0.174) GeV. The 
parameters appear in the tt distribution amplitudes are listed below. We use the values at the scale fi = 1 GeV in 
our calculation under the consideration that the heavy quark behaves almost as a spectator of the decay processes in 
our discussion at the leading order of HQET. 



02 m W3 774 uji hoo 
0.25 0.015 -1.5 10 0.2 -3.33 



Woo aio wio hoi hiQ 
-3.33 5.14 5.25 3.46 7.03 



We will work at the symmetry point, i.e., Ti = T2 = 2T, uq = 1/2. This comes from the consideration that every 
reliable sum rule has a working interval of the Borel parameter T within which the sum rule is insensitive to the 
variation of T. So it is reasonable to choose a common point Ti = T2 at the overlap of Ti and T2. Furthermore, 
choosing Ti — T2 will enable us to subtract the continuum contribution cleanly, while the asymmetric choice will lead 
to the very difficult continuum substraction fl^. 

From the convergence requirement of the operator product expansion and the requirement that the pole contribution 
is larger than 40%, we get the working interval of the Borel parameter T. The resulting sum rules are plotted with 
LJc = 3.2, 3.4, 3.6 GeV in Figs. [TJIH The numerical values of these coupling constants are collected in Table H] 



III. SUM RULES FOR THE p, uj COUPLING CONSTANTS 



The sum rules for the p (w) coupling constants can be obtained using the same apporach. Now the tensor structure 
of the decay amplitudes are a little more complicated than those of the tt meson, due to the spin of the final p {uj) 
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FIG. 1: The sum rule for 5M2Si7r/+,i/2/-,3/2 
with ujc = 3.2, 3.4, 3.6 GeV and the working 
interval 3.0 < T < 4.0 GeV. 
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FIG. 2: The sum rule for ffi?^Ti7r/+,3/2/-,3/2 
with ojc = 3.2, 3.4, 3.6 GeV. There is no working 
interval for the Borel parameter T. 




FIG. 3: The sum rule for 9M^T,7rf+,3/2.f-,3/2 
with Uc = 3.2, 3.4, 3.6 GeV and the working 
interval 3.0 < T < 4.0 GeV. 
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FIG. 4: The sum rule for .^m^Mitt/- 3/2 ^^^^ 
cOc = 3.2, 3.4, 3.6 GeV and the working interval 
3.0 < T < 4.0 GeV. 



meson. For example, the decay amplitude A4{M2 Hi+ p) can now be written as 

M{M2^Hi+p) = 2Ir) 



^Oiie qv J^a2 I „Qia2c.€ e qv 
-f- — t/j fc 



+2/??aic 



where e* is the polarization vector of the final p meson. 



9t 



It 



/2 

9M2H1P ■ 



(14) 



pi* ci2 
Hm^Hxtt yM2Sin 



sO* 



„d2 pi* /3 

i/MiTiTT yM2Mi7r !/M2Mi7r 



5e 1.13 -0.68 -0.40 -1.09 0.04 -0.46 
0.78 ~ 1.65 -1.07 0.47 -0.83 0.15 -1.67 0.75 0.03 ~ 0.05 -0.62 0.34 



TABLE I: The tt coupling constants in units of [GeV] with j the orbital angular momentum of the final pion. cjcS 
correspond to the central values of the overlap amplitudes and gs with Wc = 3.4 GeV and T = 3.5 GeV where 
9m2Hitt = 5M2ffi7i-/-,i/-,f 6tc. The ranges of gs are determined according to the uncertainty of the overlap 
amplitudes and the ranges of gs with 3.0 < T < 4.0 GeV and 3.2 < ojc < 3.6 GeV. We use asterisk to indicate the 
coupling constants from sum rules without a stable working interval. 



We define other p coupling constants as 
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(15) 



Because of heavy quark symmetry, there are only six independent coupling constants involving these decay modes. 
Their sum rules read 
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+/Jm2 [i6(w/i(|'*))[^1(mo) - 16/i[|^)['1(mo) + ^t(wo - (wAt)(uo)) - 324''(«o) - IQT^-^^^^m) 
-32T|-^'°l(no) - 16Ti-^'°l(no) - 1QtI-^'%o) + Z2S^-^^yu<,)\ ^ + /,m^ [4AP1(mo) 
-4^ti](up) + 4(mA)[i1(uo) - 32C7[3](wo) + 32(wC)[^l(wo) + 64C7W(uo) + 32VI-^'°1(mo) 
+64*1-2.01 (mo) + 64$[-2.o1(mo) + 64$[-2.o1(mo)] ^ 



(20) 



(V2 £ £ 

9M2Sipj+,y-,i^ 



'^-,3/2+^-F,l/2 
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+.fjml [{uAt){uo) - At(uo) + 16ri-i'°lK) + 16ri-^'°i(uo)] ; 
The other p couphng constants are related to the above ones by the fohowing relations: 



pi 

dAhHip — 


\/6 pi 

4 9 Ml Hop - ~ 


a/6 pi 

2 9miHip 


p2 

9m2Hxp ~ 


VQ p2 

g 9MiHip - - 


1 p2 

29M2H0P ' 


/2 

9m2Hip ~ 


g 9 Ml Hip - " 


29m2Hop ' 


9M2S1P — 


-^a^i - 

^ 9MiSoP - 


2 9miSip 


9M2S1P — 


di _ 

^ 9MiSoP - 


V6 di 

2 ilMiSip I 


d2 

9M2S1P — 


V^ d2 
-^9MiSip - 


29M2S0P ■ 



(21) 



(22) 

In our numerical analysis, the parameters that appear in the distribution amplitudes of the p meson take the values 
fromRef. [11. 



/p[MeV] /J[MeV] aJj Clp ^Ip ^Ip ^ip Q Ci U 

216(3) 165(9) 0.15(7) 0.14(6) 0.030(10) -0.09(3) 0.15(5) 0.55(25) 0.07(3) -0.03(1) -0.03(5) -0.08(5) 



9m2Hip 



p2. 
y A/2 Hip 



/2 

9M2H1P 



si. dl 

9M2S1P 9M2S1P 



d2 

9M2S1P 



Qc -0.10 -0.26 -0.62 0.16 0.32 0.30 

g -0.20 ~ -0.04 -0.40 ~ -0.16 -0.93 - -0.41 0.09 ~ 0.27 0.21 ~ 0.53 0.21 ~ 0.47 



pi* 
9M2H1U 



p2* 
9M2H1U 



/2 
9M2H1U 



si* dl d2 

yM2SiLj 9M2S1UJ 9M2S1L, 



gc -0.08 -0.24 -0.56 0.14 0.32 0.27 

g -0.13 ~ -0.04 -0.34 0.16 -0.79 0.41 0.08 ~ 0.25 0.21 ~ 0.53 0.19 ~ 0.43 



TABLE II: The p/uj coupling constants in units of [GeV] with j the orbital angular momentum of the final p/u 
meson. gcS correspond to the central values of the overlap amplitudes and 5s with Wc = 3.4 GeV and T — 2.75 GeV. 
The ranges of gs are determined according to the uncertainty of the overlap amplitudes and the ranges of gs with 

2.5 < T < 3.0 GeV and 3.2 < Wc < 3.6 GeV. 

The resulting sum rules are plotted with ujc = 3.2, 3.4, 3.6 GeV in Figs. [SJSl The numerical values of these coupling 
constants are collected in Table |lTl 

Replacing the p meson parameters by those for the uj meson, one obtains the uj meson couplings with the heavy 
mesons, which are also listed in Table HIl Here we take the following values for the parameters /(^, /J [17], and m;^: 
= 0.195 GeV, /J = 0.145 GeV, and = 0.78 GeV. 



IV. SUM RULES FOR THE 7 COUPLING CONSTANTS 



As an example, we consider the decay Mi — iJi + 7. The decay amplitude should now be expressed according to 
the total angular momentum of the final photon, namely ml, e2 • • • : 

MiMi ^Hi+j) = [(ry • e*)(e* ' Qt) - {v ' qt){e* ■ e*,)] g^l^,^ 
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(Uj=3.2GeV 
(Uc=3.4GeV 
(j,=3.6GeV 




2.5 3.0 
TIGeV] 



FIG. 5: The sum rule for Dj^f^jj^pf— ^1/2 f— ,3/2 
with cOc = 3.2, 3.4, 3.6 GeV and the working 
interval 2.5 < T < 3.0 GeV. 



FIG. 6: The sum rule for gM2Sipf+,i/2f-,3/2 with 
coc = 3.2, 3.4, 3.6 GeV. There is no working 
interval for the Borel parameter T. 





— <Uj=3.2GeV 




<^j=3.4GeV 




<Uj=3.6GeV 







2.5 3.0 
T[GeV] 



— (u.=3.2GeV 

(u,=3.4GeV 

(/)^=3.6GeV 



2.5 3.0 
T[GeV] 



FIG. 7: The sum rule for g%s^pf+,inf-,m ^"^^ 

= 3. 2, 3. 4, 3. 6 GeV and the working interval 
2.5 < T < 3.0 GeV. 



FIG. 8: The sum rule for 5m2Sip/+,i/2/-,3/2 with 
(jOc = 3.2, 3.4, 3.6 GeV and the working interval 
2.5 < T < 3.0 GeV. 



(23) 



where e is the charge of the proton. To extract 5Mi/fi7 ^^'^ SmiHi-j^ ^® consider the following correlator: 



/ 



d^xe-^'=-(7(a)|T{Jf_ ,(0)JT"_ 3(a;)}|0) = ei 

+ei 



(^t Qt - <lt 



It) 



It 



G'j^^H,-,{uj,u>'). (24) 



In the limit niQ 00, two diagrams contribute to the above correlator according to the way how the photon couples 
to the light quark. First, the photon couples to the light quark perturbatively through the standard QED interaction. 
For this diagram, we need the quark propagator in an external electromagnetic field: 

{0\T{q''{x)q''{0)}\0)F,^ = ^ du{2{l - 2u)x^j. + ie^^pal^Yx^jF'^'' {ux) , (25) 

where we have adopted the Fock-Schwinger gauge x^Ap{x) = to express the electromagnetic vector potential in 
terms of the gauge invariant -F^,,/. In our calculation, we take the u quark as the light quark involved and therefore 
Cq = I . The second diagram involves the nonperturbative interaction of the photon with the light quark in terms of 
the photon light-cone distribution amplitudes. 

Now the problem can be tackled using the same approach as in Section II and III. We list the final sum rules below: 



3/2+^-, 1/2 
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1 n_ 



T^h{^) - ^x{m) K(^o) - W^)'(uo)] tVi(7#) - 4/37 r (^0) + 2r{u,) 



+8(«/i^)W(uo) +325[0'01(uo) + 16Tf •°l(wo) - 16Ti°'°l(«o) 

^-,3/2+-^-, 1/2 



(26) 



1 



8\/6 



+V[-^'°l(wo)] ^ - 6(9(7) [^™(wo) - M)W(wo) + 8/iW(«o) - 8/iPl(wo) + ^{uh^f\u^) - 16Ti-''°](uo) 
-16Ti-^'°l(.o)] ^ 

Similarly, the electromagnetic coupling constants between doublets M and S can be defined as 



(27) 



7W(Mi ^51+7) 



el 



+e 



2£^^*«"(e* • gt) "^t^ - (e* • g*) 



m2 



(28) 



and the sum rules for g^^Si-y ^"^^ 9M%i-y read 



'2 '2 

1 J — 



7^Vi(y) - 4X(99) [-/-tI^o) - (w.^7)(^o)] Tfo{^) + 2/3^ [4^ («o) + ^r^^^uo) 



16^6 L TT 

-(wV;")(«o) +^[°'°k«o) + 8V[°'°1(mo)1 + (gg) [a(uo) - {uA){uo) + 16h^^^iuo) + 32S^-'-'Huo) 
-16T[-''°\uo) + 16Tt''°\uo) 

^-,3/2+-^+, 1/2 



(29) 



= -^e,!^^T'h{'^)+4x{qq) IM^o) - {uct>,){uo)]Tfo{'^) + 2^, [4A^°'°^{uo) - V"(«o) + (uV'")(^o) 
-(gg) \a{uo) - {uA){uo) + 16T[-^'°\uo) + 16Tt^'''\uo) 



(30) 



The 7 coupling constants of the other channels are defined as 



M{M2 Hi+j) = 2er]o 



g.aie gUg*a2 _ g.aie*e*t> 0:2 1 _ gOiia^ ^i' e* qv 



.aie qv *a2 _i_ ai£ e v 



+e??aia2|at~ 

A^(Mi ^ 5o + 7) = ei [(77 • qt){e* ■ qt) - (v ■ e*)(^] gi^^s^i ' 
M{M2 ^ 5o +7) = 2e77„,„,£«^^*«''gr53S'so7 ' 



2^a.. «-g^"=(e*.g,) 5^^^^^ 



7W(M2 ^5*1+7) 



= 2ei7?Q 



(e*-gt)- 



*ai *a2 _aia2/',* 



2 1 el 

9t rS'M2Si7 
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1.5 2.0 



2.5 3.0 
T[GeV| 



FIG. 9: The sum rule for gT/,H,J-,i/2f-,3/2 
with ujc — 3.2, 3.4, 3.6 GeV and the workmg 
interval 2.5 < T < 3.5 GeV. 



FIG. 10: The sum rule for gil^^^/-,i/2/-,3/2 
with ojc = 3.2, 3.4, 3.6 GeV and the working 
interval 2.0 < T < 2.4 GeV. 
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FIG. 11: The sum rule for gl\st-if+M2f~,3/2 
with LOc = 3.2, 3.4, 3.6 GeV and the working 
interval 2.5 < T < 3.5 GeV. 



FIG. 12: The sum rule for gTiis,jf+,i/2f-^3/2 
with uJc = 3.2, 3.4, 3.6 GeV and the working 
interval 2.5 < T < 3.5 GeV. 



+2ezr7„,„J 2 [er^q^ie* ■ qt) - q^'q^ie* ■ e?)] + [er^g^ " 9?'"He* ■ qt)] (e* ■ q,) 



9r"'{e*-e;)]q^\g^jis,^ 



(31) 



There arc simple relations among the coupling constants of a certain decay type (el, ml • • • ) within two doublets: 



ml 



9 Ml Hif 



9 Ml S 1-1 



m.2 
9MiSij 



1 



ml 
9MiHai 



2' 

V6 ,2 
'—9M2H0-1 



V6„ml 



1 



3MiSo7 



r, 9M2S0I — vOffAfaSii 



(32) 



The resulting sum rules are plotted with uJc = 3.2, 3.4, 3.6 GeV in Figs. l9lfT2l The numerical values of these coupling 
constants are collected in Table IIIII 

The parameters in the distribution amplitudes are = -(0.0039 ± 0.0020) GeV^ w^Jf = 3.8±1.8, uj^ = -2.1 ±1.0 
[Hi; K = 0.2, Ci = 0.4, C2 = 0.3, (^2 = = Ci^ = = (^ = 1 GeV) |2(|. The quark condensate at the same 
renormalization scale is (qq) = (-0.245 GeYf{q ^u,d) We adopt the value x = (3.15 ±0.3) GeV"^ fill for the 
magnetic susceptibility of the quark condensate. 
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ml 



yMiHi- 



m2 
9 Ml Si 



7 

GeV^^l 



[ GeV-i] [ GeV"^] [ GeV"^ 
Qc -0.05 0.37 0.06 0.16 

g -0.07 0.03 0.26 ~ 0.53 0.04 ~ 0.10 0.10 ~ 0.28 



TABLE III: The 7 coupling constants, g^a correspond to the central values of the overlap amphtudes and gs with 
Ua = 3.4 GeV and T = 2.75 GeV(2.2 GeV for gfj^ 

Hi-y)- The ranges of gs are determined according to the uncertaintj 
of the overlap amplitudes and the ranges of gs with T in their working intervals and 3.2 < uic < 3.6 GeV. 



V. DECAY WIDTHS 



It is straightforward to calculate the decay widths of the ( 1 ^ , 2 ~ ) I?- wave heavy mesons with the coupling constants 
extracted in the previous sections. For completeness, we make rough estimates of the decay widths of the single and 
double-pion channels with the given coupling constants, although some of the coupling constants are extracted from 
sum sules without a stable working interval, as shown in Section. Illland lTTTl The formulas for the single-pion and the 
radiative decay widths are collected in Appendix [Cj 

To get the numerical values of these decay widths, we need the mass parameters of the heavy mesons concerned. 
They are collected in Table IIVI 







0" 


1" 


0+ 


1+ 


1+ 


2+ 


1" 


2~ 




[GeV] 


1.87 


2.01 


2.40 


2.43 


2.42 


2.46 


2.82 


2.83 


ms 


[GeV] 


5.28 


5.33 


5.70 


5.73 


5.72 


5.75 


6.10 


6.11 


mo^ 


[GeV] 


1.97 


2.11 














TUB, 


[GeV] 


5.37 


5.41 















TABLE IV: Heavy meson masses used in our calcultion. We take the values of m^i, mo', i^d^, ^d[: ttt-Di, 1^1021 
"m-B, TUB*, iTiB'^, iTiBi, ^Dsj Slid tub^ from Ref. (2]| . Masses of the other heavy mesons are from the quark 

model prediction made in Ref. [22!]. 



The resulting decay widths are presented in Table [V] Here we also make rough estimates of the decay widths 
of D + rj/K^ channels assuming the flavor SU(3) symmetry. The column F^i^^j+p sums over the widths of all 
pseudoscalar meson decay channels, including the widths of Dtt'^ channels which are about half of the widths of the 
corresponding Dtt'^ channels. 

The Feynman diagram of the double-pion decays is shown in Fig. 1131 We will not present the expressions for 




FIG. 13: The double-pion decays of ZJ-wave heavy mesons via a virtual p meson. 



the widths of these double-pion decay channels due to the length of these expressions. Their numerical values are 
collected in Table IVTl together with those of the radiative decay channels. We take /pttt = 6.14 in our calculation. 
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D7r+ 








i3i7r+ 


D27r+ 




DsK+ DIK+ 






3.7 


1.3 




0.02 


2.1 


0.04 


0.3 0.1 


1.6 0.4 


13.1 




1.8 - 7.9 


0.6 - 2.9 




0.01 - 0.04 


0.3 - 8.9 


0.02 - 0.1 


0.2 - 0.7 0.1 - 0.2 0.8 - 3.4 0.2 - 0.8 


5.4 - 34.9 






4.4 


0.02 


0.02 


0.04 


7.8 


0.3 


1.2 


19.9 






2.1 - 9.3 


0.01 - 0.05 


0.01 - 0.06 


0.02 - 1.0 


1.1 - 33.5 


0.1-0.6 


0.6-2.6 


5.6 - 69.1 




B7r+ 


B*TT + 






Bnv+ 


B2'K+ 


Br) B'q 


BsK+ BIK+ 




Bt' ^ 


4.4 


1.9 




0.02 


2.3 


0.05 


0.3 0.1 


1.3 0.5 


15.2 




2.1 - 9.3 


0.9-3.9 




0.01 - 0.04 


0.4-9.6 


0.02 - 0.11 


0.1-0.6 0.1-0.2 


0.6 - 2.8 0.2 - 1.0 


6.1 - 39.0 






5.8 


0.02 


0.02 


0.05 


8.5 


0.3 


1.6 


23.5 






2.8 - 12.4 0.01 - 0.06 


0.01 - 0.06 0.02 - 0.11 


1.2 - 36.3 


0.2-0.7 


0.7-3.3 


7.0 - 77.4 



TABLE V: The widths of the single-pion channels in units of MeV, including their ranges and the values 

corresponding to the central values of gs. 







D*7r+n° 


i3o*7r+7r° 




D-y 


D*7 


D*o7 


D[j 


Dl' 


44.5 


642.3 


3.0 


1.2 


8.0 


12.2 


0.3 


0.8 




7.1 - 178.2 


241.3 - 1529.8 


1.0-8.5 


0.4 - 3.4 


2.9-15.6 


5.6 - 24.7 


0.1-0.7 


0.3 - 2.3 


^2 ^ 


1155.0 


448.3 


0.06 


3.5 


9.0 


15.7 


0.4 


0.8 




440.6 - 2728.3 


165.1 - 1092.8 


0.03 - 0.1 


1.1 - 10.0 


4.4-18.5 


6.5-31.3 


0.2 - 1.3 


0.3 - 2.3 




Bn+7T° 








S7 


B*7 


So7 


Bh 




19.1 


579.7 


2.7 


1.1 


9.4 


18.8 


0.3 


0.8 




3.1 - 76.5 


217.3 - 1382.1 


0.8 - 7.6 


0.3 - 3.0 


3.4-18.5 


8.8 - 38.3 


0.1-0.8 


0.3 - 2.4 


B2* ^ 


521.7 


414.2 


0.05 


3.1 


9.5 


22.2 


0.5 


0.8 




198.1 - 1234.1 


152.5 - 1008.4 0.02 - 0.12 


1.0-8.7 


4.7-19.6 


9.3 - 44.5 


0.2 - 1.4 


0.3 - 2.4 



TABLE VL The widths of the double-pion and radiative decay channels in units of keV. 



VI. CONCLUSION 

We have calculated the tt, p, ui, and 7 couplings with heavy mesons at the leading order of HQET within the 
framework of LCQSR. Most of the sum rules obtained are stable with the variations of the Borcl parameter T and the 
continuum threshold ujc- For the other sum rules, we can not find a stable working interval of T. The extracted vector 
meson heavy meson coupling constants may be helpful in the study of the interaction between two B{D) mesons. 

Some possible sources of the errors in our calculation include the inherent inaccuracy of LCQSR: the omission of 
the higher twist terms in the OPE near the light-cone, the variation of the coupling constant with the continuum 
threshold lUc and the Borcl parameter T in the working interval, the omission of the higher conformal partial waves in 
the light-cone distribution amplitudes of the involved mesons or photon, and the uncertainty in the parameters that 
appear in these light-cone distribution amplitudes. The uncertainty in f's and A's is another source of errors. It's also 
understood that the 1 / mg correction may turn out to be quite large for the charm mesons while such a correction is 
under control for the bottom ssytem. 

With these extracted coupling constants, we have made a very rough estimate of the decay widths of the (1~, 2~) 
Z)-wave heavy mesons. Their dominant decay modes arc M ^ T + n and Af — > iJ + 7r. The former mode is of S'-wave 
while the latter decay occurs in the P-wave but with a larger phase space. The total width of the D-wave heavy 
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mesons is roughly several tens of MeV. Therefore the (1 ,2 ) _D-wave heavy mesons are not expected to be extremely 
broad, which is helpful in the future experimental search. 
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Appendix A: The definitions of the n, p meson and the photon light-cone distribution amplitudes 

The 2-particle distribution amplitudes of the tt meson are defined as (Tsj 



Acknowledgments 



{0\u{zh,,j5d{~z)\TT-{P)) 




— z, 



1 



I 



where ^ = 2u — 1, is the leading twist-2 distribution amplitude, 4>(p^cr) ^re of twist-3. All the above distribution 
amplitudes = {(/)7r, 0p, <pa,g7T} are normalized to unity: du (j){u) = 1. 
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There is one 3-particle distribution amplitudes of twist-3, defined as [l 



\u{z)a^^j5gGaf3{vz)d{-z)\TT (P)) = 



(paPi^gi^p - PcPi^g^p - ppPt^gL + PuPfg^t,) pz) , (A2) 



(A3) 



(A4) 



"1,1 + md 

where we used the following notation for the integral defining the 3-particle distribution amplitude: 

Here a is the set of three momentum fractions a^, a„, and ag. The integration measure is 

j Da = j daddaudagS{l — au — a.d — ctg) . 
The 3-particle distribution amplitudes of twist-4 are 

{Q\u{z)-i^,-^^gGap{vz)d{-z)\n~{P)) = Pf,{paZ^ ~ ppZa) — f-^mlA\\{v,pz) + [ppgi^ ~ Pagp^)fn'mlA±{v,pz) , 

pz 

{0\u{z)'^f,igGa(3{vz)d{-z)\TT^{P)) = p^{paZp ~ ppz^) — f-^mlV\\{v,pz) + {ppg^^ - pagj^)f-^rnlV±{v,pz) , (A5) 

pz 

where Gap is the dual field Gap = \£ap-i5G'^^ . 

The definitions of the distribution amplitudes of the p meson used in the text read as [13, [3 



l{z)-i,4{-z)\p-{P,\)) = /p 



= (A) 



p^ / due^^P-^v\\{u,p^)+e% / du^^^-'g^l\u,p^) 



P-z Jo 
1 eW.- 



z,,^ / due'^P'''g3{u,fi'^) 





(0|n(z)7p75rf(-^)|p-(P,A)) = lfpm,e/'^"e^llpazp [' due^^P'^ g^^\u, fi^) , 



{0\uiz)a^,d{-z)\p-iP,X)) = if, 



Jo 

+ {Pf_,Zi, ~ pi,Zfj^)- —ml / dwe'^P'^/i /(u,^^) 

[P ■ z) Jq 

2 ^ P • -2 Jo 



(0|u(z)d(-z)|p-(P,A)) = -tfj{e^^^z)ml I due'^P-'h\'\u, n^) ^ 



(A6) 



The distribution amplitude (^|| and are of twist-2, g^j^\ 9± \ ^||*'' a-nd are twist-3 and 33, are twist-4. All 
functions (j) — {ip\\^ip±,g^^\g'f^\h^^^\h^^\g3,h3} are normalized to satisfy ^}d,u(j){u) — 1. 
The 3-particle distribution amplitudes of the p meson are defined as [TtI . [l 



'\u{z)gG^,,-1albd{-z)\p (P, A)) = /pmpPQ[pi.e^^j - p,,e\_^j]^(u,pz) 

e(A) . z 

+ /p^p [Pt^giv -pygit,]'^{v.pz) 

+ /p"-^p Pa[p^iZu-p^Z^]^{v,pz) 

{0\u{z)gGf,J-fad{-z)\p^{P)) = fpmppa[p„e''ll-p^e^llv{v,pz) 



+ /p"^p {Pi^gav - Py9ay\^{^-.PZ) 



pz 
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g(A) . ^ 

+f^^l[Pa(^''ll9fy - Ppe^llaL - Pae^!^^^^ 
rT^2 



+ — — -[PaP^,e.^llzu - PpPt^e^^^z^ - PaPue^ljiZ^, + ppp^e^l^z^^T^*' {v,pz) 
pz 



+ — — -[PaPt^e'^^lzp -Pi3Pt,e!'llza - PaPve'llzfj + pi3p^e'l'^Za]Tl^' {v,pz) , 
pz 



JA), 
JA), 



{0\u{z)gG^^{vz)d{-z)\p (P, A)) = i/jTO^[e^^^p^ - e^^^p^]S'(t;,pz) , 
(0|72(z)zgG^,(i;z)75d(-z)|p-(P, A)) = ifjmyiy - e^llp,]Siv,pz) . (AT) 

The distribution amplitudes A, V and T are of twist-3 and the other 3-particle distribution ampUtudes are of twist-4. 
The definitions of the 2-particle distribution amphtudes of the photon read asp^ 



{0\q{z)aai3q{-z)\j^^\q)) = « e x (g?) (g/jei"^^ " '7a4'^^) J due'^'^'' (l)y{u, p.) 



1 



1 /"^ 



(A8) 



(0|g-(z)7^g(-z)|7(^)(9)) - 

{0\q{zh^.l5q{-z)\j^^\q)) = 
The 3-particle distribution amplitudes of the photon are defined as [l^ 

{0\q{z)gGf,^{vz)-fa-f5q{-z)\-f'^^\q)) = Cg e fs^^ qa[q„e^ll - qf,e^ll]A{v,qz) , 
{0\q{z)gGf,^{vz)i"/aq{-z)\"f'-^\q)) = CgC fs^ qa[qi,e^ll - 9^e^^]V(i;, qz) , 

{0\q{z)gG^^{vz)q{-z)\"f'-^\q)) = ieqe{qq)[qt,e'll- q^e^ll]S{v,qz) , 
{Q\q{z)gG ^u{vz)i-i^q{-z)\-i'^^\q)) = ie^ e {qq)[qi,e'^lj^ - q^e^ll]S{v,qz) , 

{Q\q{z)aafigG^^{vz)q{-z)\-i'^^\q)) = e {qq)[qae^llg^^ - qpe^llg^^ - qae~ll9pf, + 9,3ei^^5a^]Ti(i;, qz) 

+eg e (OT)[g^e^^2.g^^ - g^e^^^.g^^ - q^e^^^gpf, + g^e^^^5^^]T2(^;, gz) 

, eqe{qq) (a) (a) (A) , (A) .rp , . 

H [qaqtie^l^Zy - qpq^^eY^Zy - qaq,,e\^Zf, + qpq^e)_^Zfj,\T3{v,qz) 

qz 

, eqe{qq) (a) (a) (a) , (A) i^ , 

H WagMe±j^/3 -qpqiiey^Za - qaqusy^zfj + qpq^e)_^Za\T4{v,qz) . 

qz 

(A9) 

Here (j)^ is the leading twist-2 distribution amplitude, A, and V are of twist-3, /i^, S*, S*, and T^j^ 2,3,4) &re 

of twist-4. 



Appendix B: The definitions of and T^' 



The functions and jr[™."] used in the text are defined as 

J"'"l(uo) = / • • • / / T{xi)dxidx2 ■ ■ ■ dxn , 
Jo Jo Jo 
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T{1 — a2 — as, a2, as) 



'^UQ—a2 



Jo J Un 



{uq - a2)J^{l - a2 - as, a2, as) 



daida2 , 



"0—0:2 

J'(l - uo, a2, uo - 0:2) , Z"^""" J'Cl - uo - "3, uo, as) 

aa2 — / das , 

Uo - Q!2 Jo "3 



J"(l - uo,a2,uo - 02) 



da2 - 



Uo — 0.2 



J"(l — a2 — Q;3,a2,a3) 



Jt*o— "2 



as 



J"(0, Mo, 1 - Uo) J"(l - Uo - Q!3, Uq, as) 



1 - Uo 



as 



a3=0 



+ /"° ^^^.^ 5[-^(l - ^2 - a3,a2,a3)/a3] 



■I 



do-s 



"3="0-«2 



^ "° das ^^'^^"'" ~ ~ Q^3iQ^2, as)] 7^0:2 



/O "3 

T{\ - Uo - a3,uo,a3) 



as 



03=0 •'o 



a2=Mo 

a[7'(l - a2 - as, a2, a3)/a3] 

aa2 ^ 

das 



a3=«o— "2 



"° /"(l - Uq, a2, Uo - a2) _ "° J'(l - Uo - as,uo,as) 

io 



+ / rfa2 
'0 



(uo - a2)2 
9[J"(1 - a2 - a3,a2,as)/3a2] 



as 



+ 



d\T{\ - a2- as, a2, as)/as] 



9as 



+ 



9[7"(0, l-as,as)/as] 



a2=Mo>a3=0 ^"3 

d{T{uo - a2, a2, 1 - uq)/ (1 - uq)] 



a2=wo,o;3=0 



+ - 



03 = 1 — Mo 



Sao 



+ / c?a2 
Jo 



d'^[J^{l - a2 - as,a2,as)/as] 



9a2 



Q:3="0— 0:2 



^^^ a^[.F(l-a2-as,a2,a3)]/aa 



as 



^[^■(l - a2 - as,a2,as)/9a2] 



as 



+ 



02=1*0 

^[^■(l - a2 - as,a2,as)/as] 



das 



+ 



T{1 — a2 — as, a2, as) 



a2=M0iO;3=0 ^^"-1^3 

T{\. — a2 — as,a2,as) 



a2=W0,O!3=0 



+- 



0^2 — 'Uojag— 



r"o - a2 - a3,a2,a3)/a3] 

+ ^ '^"^ M 

+ ^^,^ ^[-^(1- "2 - a3, a2, a3)/ai] 



7 

7' 

JO 

■ / das 

Jo 



a% 



a3="o — "2 



a2=Mo>03 = l-«o 



9as 



"3="0— «2 



^ "° 9J^(1 - a2 - as,a2,a3)/(?a2 



/ / J^(l — a2 — as, a2, a3)rfa3rfa2 
io Jo 

^ n f'-"' (up - a2).F (l - a2 - as, a2, as) 

•JO -J -Uo — 0:2 



as 



dasda2 ■ 



2 7) io 

r-o /■i-"2 _ c,2)2^(l - a2 - as, a2, as) 



J^(l — a2 — as, a2, as)rfasrfa2 
- / / — — " — — "''""'^'""'' rfasda2 , 

^ Jo Jun 
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J^I~^'°l('Uo) =11 / — a2 — x,a2,x)dxdazda2 

Jo Jo Jo 

+- / / a3iF{l — a2 — az,a2,cx.3)da3da2 

^ Jo Jo 



1 



«o /•!— a;2 



(uo - a2y 



J^{1 — Q!2 — 0^3, 0:2, o:3)da3da2 , 



+ 



2 Jo Juo-a2 "^3 
Mo rtio—a-2 pas 

I II T{1 — a2 — Xi,a2,xi)dxidx2dazda2 

Jo Jo Jo 

^ Jo JO Jo 

^ /■«o /•"o-a2 
6 

{uo - 02)' 



xT{l — a2 — X, a2,x)dxda3da2 



Jo 

110 ;'l-a2 „ \Z 



as 



— Q!2 — Q!3, a2, 0L-^da'ida2 ■ 



Appendix C: The formulsis for the single-pion and the radiative decay widths 

The single-pion decay widths can be expressed in terms of the 3-momentum q of the final pion: 

_ 1 ruHo^f pi ^2l |3 

1 1T^Hi ( ^pl n2|_|3 



r(Mi Ho'k"' 
r(Mi i^l7r^ 

r(Mi ^ SiTT^ 

r(Mi -> Tl7^^ 



r(Mi 
r(M2 



,+ 



r(M2 ^ 5o7r^ 

r(M2 SlTT^ 



r(M2 
r(M2 



■Tl7r^ 
• T27r^ 



12^ mM/^"^^^^"^ ' 



= d2 X2| |5 . 1 1 

- Stt mM, ^^^^^^^^ ^ 24^ mM.mT, ^^^^^^^^ IStt mM^rriT, 

1 Ll!!ZLc„d2 )2| |5 I 1 ^ / d2 -,21 |7 

247rmM2^^^=^^"^ ' ' mM^mn}^''-"'^' ' 



^2|„|5 



(Bl) 



rr+^ - ^ "^■S'o (-„d2 n2| |5 

- "^-gi f„rf2 -.21 |5 

- J_Z!!Zi_('fl'^2 -,2, |5 

- in \9M2TnT) |y| ) 
""''^ ~ ^ (5M2T27r)^k| + ■^ZZ; Z; (SMsTstt)^!?!^ " 7TZ " Z;3~ (5M2T2 7r)^ kl 

zTT mM2 <j7r mM2 '^T2 4o7r mM2 '^t2 

, 7 9m2T27t9m2T2-k i I 7 "^T2 ^ „d2 ^2|„|5 , 2 9M2T27r9M2T2^ |„|7 

+ 307r mM2mT. ^ SOtt mM2 ^'^'''^^^"^ 457r mM2m3^ 

The above formulas for the decays M(i^2) — ^ -^(0,1) + ^'^'^ also valid for the decays M(i_2) ^ -^(0,1) + 
except for an extra 1/6 isospin factor in the case of r] decays. 



Similarly, the radiative decay widths read as 



2 2 1 2 ?n 1 c2 
L{Mi^Hi'y) = [dMiHi-y) 111 + T7i [gMiHi-y) 9 + 



6tt mMimsi 127rmMims;^ 

rIM ^ M ^\ - ^^'"1 ^2|_|3 , 1 ^„ml \2|„|5 , ^M^ji.-'Jn^nn I 

i(M2-)>iii7j - — — — (5M2ffi7j ^ +Tn~::;; \9m2Hij) m +77; — — — 

+ ^^7; ^(5'M2Hi7)^|9r + 77; ~; {9M2Hi-ifWt ' 

207rmM2 iOTTTTlMi^Hi 

2 /7"^2 ^el 2 1 
e t/M2Si7^M2Si7 |„|7 I e f„m2 ^2l |7 

"■5 9 + Tri VffM2Si7j 9 • 



